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We describe a rich family of binary variables statistical mechanics models on a given planar 
graph which are equivalent to Gaussian Grassmann Graphical models (free fermions) defined on the 
same graph. Calculation of the partition function (weighted counting) for such a model is easy (of 
polynomial complexity) as reducible to evaluation of a Pfaffian of a matrix of size equal to twice the 
number of edges in the graph. In particular, this approach touches upon Holographic Algorithms 
of Valiant 1-3 and utilizes the Gauge Transformations discussed in our previous works [3HS] . 
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I. INTRODUCTION 

This paper rests on classical results derived and discussed in statistical physics and computer science. Onsager 
solution 9 of the two-dimensional Ising model on a square grid, its combinatorial interpretation by Kac, Ward [10] 
and Vdovichenko [TT], and the relation between the Ising model and the dimer model, established independently by 
Temperley, Fisher jl21 113| and Kasteleyn |14[ 115) , have aimed mainly at analysis of phase transitions in parametrically 
homogeneous infinite systems. However, the algebraic and combinatorial techniques used in these papers were intrinsi- 
cally microscopic, and thus suitable for analysis of a much broader class of problems, e.g. these that are parametrically 
inhomogeneous (glassy) and are formulated on finite planar graphs. In fact, the approach was extended |13L 116) to 
inhomogeneous versions of these models formulated on an arbitrary finite planar graph of size N. It was shown that, 
under these conditions, the partition (generating) function of an arbitrary Ising model (however without magnetic 
field) or of the dimer model are represented by Pfafhans (or determinants) of the properly defined N x N matrices. 
Note that computation of a Pfaffian/determinant requires utmost steps, while calculation of the partition function 
for a generic statistical physics model is a task of likely exponential complexity. More accurately, this is the task of 
the #P complexity class according to computer science classification, started with the classical papers of Cook [17] 
and Karp [IS . The feature of the generic statistical physics model makes the easiness of the planar Ising and 
dimer models surprising and exceptional, especially in view of the results of Barahona [19] , who showed that adding 
magnetic field (even homogeneous) to the planar Ising model immediately elevates the weighted counting problem 
(this is the partition function computation) to be of the #P-complete class, i.e. of the complexity necessary and 
sufficient for solving any other problem that belongs to the #P hierarchy ^. Similarly, adding monomers turns the 
dimer-monomer model on a planar graph (or even on a planar bi-partite graph) into a problem that belongs to the 
#P-complete class [21'^. 

One wonders if the easiness of the dimer and Ising models on planar graphs is a lucky exception or may be it is, in 
fact, just a little top of yet unexplored iceberg? A new approach of Valiant [IHSl) coined by the author "holographic 
algorithm", shed some additional light onto this question. In [P-S^^ Valiant described a list of easy planar models 
reducible to dimer models on planar graphs via a set of "gadgets" . The gadgets were of "classical" and "holographic" 
types. A classical gadget describes an elementary graphical transformation, applied to the variables defined on a 
graph element, e.g. edge or a region, that preserves a one-to-one correspondence between the configurations of the 
original and transformed models. An example of a classical gadget would be an approach used by Fisher [13 (point- 
to-triangle transformation) to map the Ising model onto a dimer model. A holographic gadget of [IH3] involves a linear 
transformation of the parametrization basis for the binary variables, so that the solution fragments of the original 
and mapped models would be in a mixed relation when certain sums, rather than individual elements of the original 
and derived formulations, are related to each other. The freedom in choosing an arbitrary nonsingular basis for the 
holographic transformation was discussed in ^I]-[3] , however, it was explored there in a somewhat limited fashion. One 
emphasis of [IH3] was on generating a polynomial time algorithm (via reduction to a determinant) for a number of 



^ The #P, pronounced sharp-P, and #P-complete classes were introduced by Valiant in [20], who studied the complexity of calculating 
the matrix permanent. 

^ Dimer and monomer-dimer models from statistical physics are called perfect matching and matching models, respectively, in computer 
science. 
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problems for which only exponential time algorithms where known before, such as counting the edge orientations on 
a planar graph of maximal degree 3 with no nodes containing all the edges directed towards or away from it. This 
model belongs to the class of "ice" problems studied earlier in statistical physics p2H24| . We will actually discuss 
this example in details later in Section |V C| A class of easy factor-function models, stated in terms of binary edge 
variables on planar graphs of degree not larger than three, was also discussed in [&. These models were stated in 
terms of a set of algebraic equations, one per vertex, constraining the factor-functions of the model. 

In this paper we present a family of statistical models, stated in terms of binary variables defined on the edges of 
a given planar graph, with the partition functions reducible to Pfafhans of square matrixes of the size equal to twice 
the number of graph edges. The general model of this class (easy on a given planar graph) can be defined in two 
consecutive steps: (a) construct an arbitrary Wick Binary Graphical Model (WBG) on the graph, and (b) further 
apply an arbitrary Gauge transformation of the type discussed in our early work devoted to the Loop Calculus (LC) 
approach Of the two steps, both required to characterize this class of easy binary models on a given graph, 

the description of the easy WBG models on a given planar graph and consecutively showing that it is Pfafhan (or 
det)-easy, is the most novel contribution of the paper. Therefore, we find it useful to informally describe a general 



WBG model already here in the Introduction. (See Section II A for the formal description.) 

A WBG model is defined in terms of painting the graph edges (each edge can be colored or not) with the "even 
coloring" requirement (for each vertex the number of colored edges adjusted to it should be even). The weight of 
an allowed configuration (an even coloring) is given by the product of the vertex weights/contributions. A vertex 
contribution is equal to unity when no edges adjusted to the vertex are colored. The vertex weights are introduced for 
all pair colorings, defined as painting of pairs of edges adjusted to the vertex. Thus, the total number of unconstrained 
parameters associated with a vertex is Ny{Ny — l)/2, where Ny is the vertex degree (valence), i.e., the number of 
edges adjusted to the vertex. The weights of all higher-degree colorings (4, 6, • • • edges of the vertex are colored) are 
not independent, but rather explicitly expressed in terms of the pair weights according to a simple rule, illustrated 
below using an example of a vertex of degree four. Assume that the edges of the degree four vertex are numbered 
clockwise 1, 2, 3, 4 with the corresponding weights of the pair-wise colorings being wi2, W13, W14, W23, W24, Then, 
the weight of the four-edge coloring is W12W34 — Wi3i/;24 + Wi4^^23- We refer to the described rule and the corresponding 
model as the Wick rule/model, since modulo signs, the expression for the higher-order weights is represented by the 
Wick's decomposition of a higher-order (even) correlation functions of a Gaussian random field in terms of the pair 
correlation functions (covariances) . The sign rule for an individual contribution (that corresponds to a partitioning 
the even number of 2k colored edges into k pairs) to the weight is defined as follows: (i) represent the edges by 
points on a circle (according to their cycling ordering), (ii) connect the points that represent the paired edges with 
lines, (iii) count the total number of the line crossings inside the circle, (iv) choose plus/minus sign if the number is 
even/odd. Note, that the Ising model, the dimer model and some other planar easy models discussed in [iHS], can 
be represented, after a gauge transformation, in terms of a special case of the general WBG model. For example, an 
WBG representation for the Ising model is directly evident from the respective high-temperature expansion/series. 



(See [25] and discussion of Section VD 
In this manuscript we show that 

• The WBG model with an arbitrary choice of the edge associated weights, w, is equivalent to corresponding 
Fermion Gaussian model, referred to as the Grassman Gaussian Graphical (G^) model, and defined on the same 
graph. With a suitable choice of the so-called complete orientations for edge and triplets of neighboring vertexes 
contributing definition of the G'^ model (and thus called in the following suitable complete orientation of the 
graph, or when it is not confusing simply -suitable orientation), the two models are equivalent in the sense 
that their partition functions are equal to each other and also equal to an explicitly defined PfafRan of a skew- 
symmetric matrix of the size twice the number of edges in the graph. Obviously, computational complexity of 
the object is polynomial in the graph size. This main result of the manuscript is formally stated in the Theorem 
|IL3 

• Application of an arbitrary and graph-local gauge transformation (described in terms of a pair of 2 x 2 matrixes 
per edge orthogonal to each other) to the WBG model generates another det- easy model. 

This paper builds upon, and in a certain sense extends, the classical results of Kasteleyn [121 ES] and Temperley, 
Fisher [T2 |T3] developed for a planar dimer model. One key ingredient of the Kastelyan approach was the con- 
struction of edge-orientations. It was shown that for a special choice of the edge-orientations, so-called Kasteleyn 
edge-orientations, the partition function of the dimer model becomes a Pfafhan of a skew symmetric matrix built 
from element-by-element product of the edge weights and the ±1 Kasteleyn edge-orientation factors. Our scheme 
also requires defining graph-local orientations, however of an extended character, where not only edge-orientations, 
but also triplet-orientations are needed. The combination of the edge-orientations and triplet-orientations defined on 
a given graph will be called complete orientations. A special choice of the complete orientation, the suitable complete 
orientation, in a sense generalizes the notion of the Kasteleyn edge-orientation. See Section In] for accurate definitions. 



3 



Another important graph-related structure utilized in this paper is the so-called even generalized loops, decomposed 
into combinations of embedded orbits on the graph. For any finite graph these graph objects are finite, which allows to 
represent relevant partition functions in terms of a sum over a finite number of contributions associated with embedded 
orbits. At this point it also appropriate to refer a related, yet different, approach that operates with (random) walks 
and associated immersed orbits, that has been formulated for the Ising model in the classical papers of Kac, Ward 
[To] and Vdovichenko [11]. ( See also [23 built upon [TOllTl] and describing how to use the Kac, Ward-Vdovichenko 
method to solve the dimer model.) In contrary with the case of the embedded orbits, which play a key role in our 
manuscript, the number of immersed orbits is infinite even on a finite graph (since the length of the immersed orbits 
is not restricted), and thus the full partition function for a finite graph (which is also equal to a finite determinant) 
is represented in terms of an infinite product (not sum!) over the equivalence classes of the walks over the cyclic 
permutations. In this context the language of the Grassman/Fermion models extensively used in our mansucript is 
more appropriate for the finite graphical objects, i.e. embedded rather than more general immersed orbits, because 
of the nilpotent feature of the Grassmann variables. 

We would also like to note that in this manuscript we did not give exhaustive description to the full class of det-easy 
problems on a given planar graph but described its rather broad, but still incomplete, subset. This incompleteness is 
in fact illustrated on a "counter-example" of the so-called # X-matching model of Valiant [IHS] , which is reducible to 
a det-easy model, however, via a set of more general (extended alphabet) gauge transformations, rather than binary 
gauge transformations discussed in our manuscript. 

The manuscript is organized as follows. Section [IT] starts with a brief description of the key graphical objects 
used through out the manuscript and then it is split into three Subsections. Sections II A| and |II B| define the Wick 



Binary Model (WBG) model and the Grassmann Gaussian Graphical (G ) model respectively, while Section II C 



ed 
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states the main result of the manuscript concerning relation between the two models under the proper, so-cal 
suitable, choice of orientations entering definition of the G^ model. The respective Theorem is proved in Section 
in a combinatorial way, while discussion of alternative topological considerations is given in Appendix |Xj Section 
is devoted to constructing efficiently a suitable orientation for a given planar graph. The construction of a suitable 
orientation is split into two steps. First step, discussed in Section [IV A[ describes a procedure of building an extended 
graph, defining a Kasteleyn orientation on it and then generating an induced orientation on the original graph. On the 
second step, described in Section [IVB[ wc show that the induced orientation is a suitable orientation, thus resulting 
overall in a simple algorithm for evaluating the partition function of any WBG model on the planar graph. Section [V| 
discusses some details and examples of the gauge transformations, which transforms a WBG model to other det-easy 



models. The Section is broken in four Subsections. Section VA based on materials from [lHH]j contains a brief 



review of the gauge transformation procedure, as well as examples that illustrate how the gauge transformations 
reduce the dimer model, ice model and Ising model on a planar graph to the general WBG models, given in Sections 
|VB|VC|VD| respectively. The aforementioned "counter-example" , which is det-easy yet requires a more general 
reduction procedure, is discussed in Section jVEj An alternative generalization of the gauge transformation approach 
is discussed in Appendix |B[ Section VI concludes the manuscript with a brief summary and a discussion of future 
challenges. 



II. KEY OBJECTS, MODELS AND STATEMENTS 



We start describing briefly some common definitions used in this Section but also through out the paper. 

Consider a finite connected non-oriented graph, Q — {Go,Gi), with Qq and Gi C {a G 2^" |Qa,rd(a) = 2} being 
the (finite) sets of graph nodes/vertices and (undirected) edges/links, respectively. For a,b G Go we will often use a 
notation a ^ bto describe the relation {a, b} £ Gi, i.e., "a and b are connected by an edge" . A directed edge (that goes 
from node a £ Go node b € Go) is an ordered pair (a, b) & Go ^ Go with {a, b} & Gi- There are exactly two directed 
edges (a, 6) and (6, a) that represent (are associated with) a non-directed edge {a,b} € Gi- triplet a — > & — > c is an 
ordered set (o, b, c) € Go ^ Go x Go with {a, b}, {b, c} € Gi and a ^ c. With a minimal abuse we will use a convenient 
notation (a, &) G G\ or even (a, &) € ^ to describe the relation {a, h\ gG\, and (a — )■ 6 — )■ c) G 5i or (a — > 6 — )• c) G C/ 
to describe the relation "a — > 6 — > c is a triplet in fj" . 

All graphs discussed in this manuscript are planar, i.e., it can be drawn on the plane in such a way that its 
edges intersect only at their endpoints. 

The graph degree/valence is the maximal degree/ valence Sa{G) among its vertices a. A graph is called homogeneous 
if all its vertices have the same valence. A connected graph G is called irreducible if after withdrawing any single 
edge the graph stays connected; it is called reducible, otherwise. Starting with a reducible graph and withdrawing 
one-by-one the edges that increase the number of connected components we arrive at a disjoint union of irreducible 
graphs, referred to as the irreducible components of G- The edges withdrawn in the above procedure are referred to 
as the connecting edges. Note that if we contract the irreducible components of G to points, the resulting graph will 



(a) (b) (c) (d) 

FIG. 1; Illustration of the Wick model construction. Four graph illustrates four representative terms (of the total of (2n— 1)!!) 
for different 6-vertex contributions into Eq. Q, correspondent to (a) W^ffi^g] j2,5],[3,4]} ~ WigW25W34, [zero crossing], (b) 
W^{[i,2i,[3,5],[4,6]} = -W12W35W46 [one crossing], (c) W{[i,3],[2,5],[4,6]} = VF13W25VK46 [two crossings], and (d) W{[i,4],[2,5],[3,6l} = 
— W14W25W3Q [three crossings], where the upper index of the 6-vertex is dropped to simplify the notations. 

become a tree, whose edges are represented by the connecting edges of G- Also note that an irreducible graph does 
not have vertices of degree one. A linear graph C of length / is a graph with the set of vertexes and edges being 
'Co = {aoi • • • 7 a/} and £1 — {{ao,ai}, . . . , {a;_i,aj}} respectively. A circular graph S of length I has the following 
set of vertexes and edges: Sq = {oi, . . . , a/} and Si = {{ai, 02}, . . . , {fli-i, a/}, {a;, ai}}. For an irreducible graph Q 
of degree three or higher a linear segment is a linear subgraph C <Z Q with the degree of vertices (in Q) ai, . . . ,ai_i 
being two, and degree of oq, o/ exceeding two. For a circular segment the vertices a2, . . . , a; have degree two, whereas 
the degree of ai exceeds two. All irreducible graphs of degree two are circular. An irreducible graph of degree three 
does not have circular segments. 

The layout of the the material in the remainder of the Section is as follows. The Wick Binary Graphical (WBG) 
model on a planar graph and the auxiliary objects, e.g. even generalized loops and intersection index, required for the 
model definition are introduced in Section [II A| Then in Section |ITB] we proceed to defining the Grassmann Gaussian 
Graphical (G^) model and respective auxiliary graph objects, the complete (combined edge- and triplet-) orientation. 
Finally, in Section fll C| we define the notion of suitable complete orientation and state the main result of the manuscript 
describing relation between the WBG and the G"^ models defined on the same graph. 

A. Wick Binary Graphical Model 

Consider a graph of an arbitrary degree and introduce the following notations 

Definition II. 1 (Generalized loop. Even Generalized loop= Z2-cycle.) Generalized loop ofQ is a subgraph of 
Q which does not contain vertices oj degree one (an empty set is also a generalized loop). Even generalized loop, or 
'Z,2-cycle, is a subgraph with all vertices of even degree. 

Remark(s): Obviously, a simple cycle of G is also a generalized loop of G- The term Z2-cycle establishes a connection 
with a language common in algebraic topology. A Z2-chain, c € Ci{G',^2), is defined by a set of Z2-numbers c = 
{cab = = 0, l|{a, 6} G Gi}- We further introduce the standard boundary operator with its value at a vertex a £ Go- 
d{c)a = J2b~a^<ib{^od 2). A chain c is called a Z2-cycle, if it has no boundary d{c) — 0, i.e., Va G Go- d{c)a — 0. 
Note that each configuration of a binary edge model is obviously represented by a Z2-chain from Ci{G\'^2)- The 
requirement for and EBG model to be even means that the allowed configurations are represented by the Z2-cycles on 
the graph, i.e., cycles which belongs to Zi[G]'^2) = {c e Ci(C/; Z2)|9(c) = 0}. Even generahzed loop, as a subgraph 
of G associated with a Z2 cycle 7, is constructed of edges with nonzero contributions, 7a6 = 1. Therefore, in what 
follows we may think of 7 as either a Z2 structure on the graph or as of respective subgraph. 

We will also need the following graph/ vertex- local notations. 

Local partitions and intersection index. For a finite set A with even number card(A) = 2k of elements let 
P{A) be the set of all partitions of A in k distinct non-ordered pairs. An element ^ € -P(^) can be viewed as a 
set of k elements, represented by pairs {a, a'} e ^ (with a, a' e A) that form the partition. We further note that, 
since our graph G is planar, each vertex b is equipped with a cyclic ordering (we choose counterclockwise) of its 
adjacent vertices (neighbors) a ^ b. For a vertex b we can place its neighbors a ^ 6 on a circle S standardly 
embedded into a plane, according to their cyclic ordering, and denoting by 1^ — I{a,c} a segment of a straight line 
that connects a and c we introduce the modulo-2 intersection index Jq, • 1^' G Z2, associated with two distinct pairs 
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a = {a, c} and a' = {a' , c'}, as the number of intersections of with Obviously we can associate with a partition 
^ e P({ai, . . . , a2k}) of an even subset of neighbors of 6 a Z2 intersection index by 

7Vf,({ai,...,a2fc};e) = • /«' (mod 2) . (1) 

We are now ready to define one of our basic models. A binary statistical model on the graph is defined in terms of 
its partition function, describing a weighted sum over all allowed configurations. 

WBG model. The partition function of the Wick Binary Graphical (WBG) model, dependent on the vector of 
edge- weights W, is given by 

Z^SGiW)^ Y: n <l...a..>.{a|a.,.„. = l}> (2) 

7={7a6}e-2i(e;Z2) beGo 

^ ^ib) ^(fc) = , ^Njv,({ai,...,a2.};«) rr (3) 

CeP({ai,...,a2fc}) ae? 

where in the first line the summation goes over all Z2-cycles, i.e. even generalized loops (contributions associated 
with all odd generalized loops are assumed zero), of Q , and subsequent product is over all vertices of nonzero degree 
contributing the even generalized loop. The second line in Eq. ([s]) represents explicit expression of the higher vertex 

weight with fc > 2 in terms of the lowest vertex weights Wa%2 : associated with edges, where the latter represents the 
case k = 1. 

Remark(s): (a) Due to symmetric nature of the weights we can write Wa ^ = ^{a"* a } ~ ^0.% — Wa'^\. for 

(b) Overall, the expressions given by Eq. ^ can be interpreted as follows. Each contribution W^''l_^ to the sum is 
determined by a partition of the set {ai, . . . a2k} into distinct pairs ai, . . . a^. The contribution is given by the product 
Wi\^ . . . Wi''^ with the overall sign in front to be plus or minus, depending whether the total number of intersections 
between the segments la^ is even or odd. To make the construction, required to define the Wick model, transparent 
we illustrate it for S^{a) / 2 = 3 in Fig.fl] 

(c) We call the model of (2]3) Wick model to emphasize relation to the Wick theorem/rules used in the quantum field 



theory to express high-order correlation function via second-moments (covariances) in the case of Gaussian statistics. 
More specifically, the (2/c)-th order contribution ^^^j to the weight function is given by the Pfaffian of a 

(2k) X (2fc) skew-symmetric matrix, with the matrix elements given by the binary weig hts W^%, that can be viewed 
as a Gaussain integral over Grassman variables. 

(d) The introduction of the set of constraints given by Eq. ([s]) is motivated by our desire to establish term-by-term 
relation between the yet to be defined G'^ model and the discrete model 

B. Orientations and the Gaussian Grassmann Graphical Model 

As in the preceding Subsection we will start introducing key objects (orientations) required to define the G'^ model. 

Definition II. 2 (Edge-orientation.) A graph edge- orientation cr is a particular choice of a directed representative 
(a, 6) for any undirected edge {a,b} € Qi, or equivalently a set of numbers Uab = ±1 associated with directed edges 
{a,h) with the skew- symmetric condition Uab = —o'ba- 

Definition II. 3 (Triplet-orientation. Complete orientation.) A graph triplet- orientation <; is a set of numbers 
^ac"* = ±1, asso ciate d with the triplets a b ^ c with zH'c — ~^ca ■ A set (<t,^), with cr being a graph edge orientation 



(see Definition II. 2) is called a complete orientation . 



^ In the following we will use the term orientation to describe complete orientation. 
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Remark(s): Note that there are no triplets a ^ 5 — c, associated with the nodes b G Qq oi valence (degree) one. 
Therefore, the triplet-orientation components c^a?, associated with the nodes b of degree one are represented by an 
empty set. 

G'^ model. The partition fmrction of the G"^ model, dependent on a complete orientation of the graph, (cr,^), 
and also dependent on the vector of edge-weights W, is defined in terms of the following integral (sum) over the 
Grassmann anti-commuting, variables tp: 

/exp U E 'fab^c^Wl^'^ipac] expnY.{a,b)eg,Vab'Jab'Pba) U diPab 

Zg3{(t,<;;W) = ^ -r- r , (4) 

/exp ( 2 E(a,6)eei '/'abCTafcy'fca) 11 difab 

(a,b) 

where V(a, 6), (c, d) G Gi fabfcd — —^cdfab- There are two independent Grassmann variables per edge (one per 
directed edge), and the Berezin/Grassman integrals in Eq. Q are defined according to standard Berezin rules, 
V(a, 5) e CJi : / dipab = 0, / (pabdifiab = 1- 

Remark(s): (a) The partition function of the G^ model can also be interpreted as a result of the following Gaussian 
statistical average (expectation value) over the Grassmann variables 

\ \ (b^a^c)eei 

(b) The denominator in Eq. Q, which is ±1 by construction, is introduced for convenience to enforce the normalization 
condition, ZQ3(^;cr;0) — 1. 

(c) Significance of the G'^ model, as of any other Gaussian Grassman model, is in the fact that its partition/generating 
function is a Pfaffian. Indeed, the denominator in Eq. ^ is the Pfaffian of the skew-symmetric matrix (of the size 
twice the number of edges in the graph) with the following elements 



bc'^bc^^ i = (a, b) k j ^ (a, c), where 6 7^ c ~ a. 



'■^ \ (Jab, i = ia,b), k j = (b,a). 

Thus calculating the partition function of the G"^ model is easy, as requiring at most 0(|t?ip) operations, 
(d) When the graph is a tree, the numerator in Eq. Q does not depend on W , and therefore the partition function 
ZQ3(cr,^; W) = 1 is trivial. Nontrivial dependence of the partition function on W appears due to the loops in the 
graph, as will be discussed in details later in the manuscript. 



C. Equivalence between the WBG model and the model 



The main subject of this Subsection (and also of the manuscript) will consist in establishing equivalence between 
the WBG model and the G'^ model for a given (arbitrary) set of edge-weights, W . However, the partition function of 
the former model depends of W only, while the partition function of the later model depends on both the edge- weight 
vector W and on a complete orientation (cr, <;) of the graph. Therefore, the aforementioned equivalence between the 
models is established only for some special choice(s) of the complete orientation which is detailed in some number of 
forthcoming definitions preceding the main theorem/statement. 

Definition II. 4 (Closed, immersed and embedded v^falk. Immersed and embedded orbit. Oriented edges 
and triplets of the walk/orbit.) A closed walk C of length 1{C) on a graph Q (not necessarily planar) is an ordered 
set C = (cq, . . . , c„) of nodes Cj G Qq with c„ — Cq. A closed walk is called immersed if Cj_i =/= c^+i, Vj = 1, . . . ,n, 
i.e., the walk does not involve backtracking events (hereafter we use a natural agreement c„+i = ci). A closed walk 
is called embedded if Cj ^ Ck, Vj ^ fc, i.e., the path does not have self-intersections. An equivalence class of closed 
walks with respect to the cyclic permutations of the path nodes is called an orbit. An immersed orbit and an embedded 
orbit is an equivalence class, represented by an immersed closed walk and an embedded closed walk, respectively. The 
ordered pairs {cj-i,Cj) and triplets (cj_i, c^-, Cj+i) for j = l,...,n are called the oriented edges and triplets of C, 
respectively. 
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Remark(s): Obviously any embedded closed walk (or orbit) is immersed. The definition (II.4) assumes that orbits 
are oriented, however we will also be discussing undirected orbits as equivalence classes of orbits with respect to 
orientation change, sometimes without additional clarification (when it does not cause a confusion). An undirected 
embedded orbit is synonymous to a simple cycle. The terms immersed and embedded are borrowed from topology and 
are used here based on the fact that immersed and embedded closed walks can be viewed as discrete counterparts 
of the loops, immersed and embedded, respectively into a plane. An immersed loop into a plane, or equivalently an 
immersion S-» of a circle into a plane is a smooth closed trajectory with everywhere non-zero velocity. This 
analogy is detailed in Appendix |X] and applied to unveil the topology that stands behind the equivalence of the WBG 
and models. 

Sign of an immersed orbit, associated with a given complete orientation. Given a complete graph 
orientation (<t,^), we associate a sign e{C) = ea-,^{C) = ±1 with any immersed orbit C by 

/ \ i(C)-l i(C)-l 

<c)^[ n n -'^"^ n 4^-1^... n -^.'=..- (6) 

\(b^a^c)eC / (a,b}£C j=0 j=0 

Remark(s): Note that the binary function e is obviously invariant with respect to cyclic permutations of C, therefore, 
the signs e(C) are defined for oriented immersed orbits. A change of orientation changes the sign of each factor in 
the products in the rhs of Eq. (|6|, and since the number of the tr-factors is the same as the number of (^-factors, 
s{C) is invariant with respect to the orientation changes. Therefore, e is actually correctly defined for a non-oriented 
immersed orbit as well. 

Decomposition of a cycle. A generic Z2-cycle can be decomposed into a set of immersed orbits in a variety 
of ways. It is easy to realize that each decomposition can be labeled by a set of variables ^ = {^ajoGeo' referred to 
as partition data where for a £ ^ denotes a partition of the set of the adjacent edges {a, 6} G 7 (or equivalently 
adjacent nodes 6 G 7) into distinct non-ordered pairs. Formally can be viewed as a set whose elements are given 
by the aforementioned distinct pairs. The partition data ^ naturally provides a set {6*1(7,^), . . . , C„(.y_|) (7, £)} of 
immersed orbits that decompose 7. The partition variables have been already used earlier to express the higher 
vertex weights of a Wick binary model in terms of the lowest vertex weights, see Eqs. ([2| and ([3|. 

Definition II. 5 (Decomposed cycle. Decomposition Components. Total self-intersection index of a 
decomposed cycle.) A pair (7,^), where ^ is a Z2-cj/c/e and ^ is the partition data, associated with its vertices 
is called a decomposed cycle. The set {Ci(7,^), . . . ,C„(..y_l)(7,4)} is called a decomposition of j associated with ^. 
The immersed orbits Cj{'j,^) are referred to as the decomposition components. The total self-intersection index of a 

decomposed cycle iV(7,^) G Z2 is defined as a sum A^(7,|) = J2aei -^ail,^) ^ ^ae-y I^{,^ °'f^° ' ^V^^^ 



vertices, where the (elementary) intersection index • la' was defined in Subsection\II A 



Remark(s): (a) The definition of the local intersection index Naij, ^) is fully consistent with Eq. ([T]), since A^a(7, ^) = 
A^o(7a,Ca), where 7a is the set of neighbors of a that belong to 7. 

(b) Note that the decomposition components are immersed orbits of a special type: they do not have self-intersecting 
and intersecting edges, and all intersections and self-intersections (if any) occur at the vertices only. In particular, the 
components generate Z2-cycles that will be denoted in the same way, and the decomposition can be also represented 

as7 = ES';''^Q(7,0- 

The just introduced notion of a decomposed cycle allows a convenient definition of a suitable orientation. 

Definition II. 6 (Suitable Complete Orientation.) A (complete) orientation is called suitable if for any decom- 
posed cycle (7, ^) the following relation holds 



(_1)9(7,€) ^ 1^ (7) 



where the binary function q(-f, ^) is defined by 



»(7.€) 

(_1).(7.C) ^ (_i)^(7.«) Yl (-£(Q(7,0)) ■ (8) 

We are now in the position to state the main theorem. 

Theorem II. 7 For any connected planar graph (of an arbitrary degree) the discrete-variable Wick model and 
fermion model are term-by-term, i.e. fully, equivalent to each other, i.e., Zq3(<;, (t;W) ~ ZiYBci^) ; for any W, 
provided (cr,^) is a suitable orientation of Q . 
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III. PROOF OF THE WBG AND G^* EQUIVALENCE 



The main task of this Section is to sketch the proof of the Theorem 11.7 We start by making the following (almost 



obvious) statement that reduces the proof of Theorem |II.7| to the case of irreducible graphs Q. 

Proposition III.l (i) The partition functions Zy^rBcO^) and ZQ3(<;,cr]W) of the WBG and models both satisfy 
the factorization property: the partition function for Q is given by the product of the partition functions of the corre- 
sponding models restricted to the irreducible components of Q. (ii) An orientation {(t,<;) on Q is suitable if and only 
if its restriction to any irreducible component of Q is suitable. 

Proof: For the WBG model both the factorization property (i) and the feature (ii) follow immediately from the fact 
that any Z2-cycle on Q is restricted to the disjoint union of its irreducible components. The factorization property 
for the G'^ model is obtained by consecutive integrations of the numerator and denominator of Eq. (4]) over the pairs 
difabdipba of Grassmann variables related to the connecting edges {a, b}. Each integration is performed by representing, 
exp((/3afc(pfca/2) = 1 + tpabfba/'^- Substitution of the first term into the numerator/denominator of Eq. Q results in 
a product of two Grassmann integrals of even functions over an odd number of integrations, which yields zero. 
Substitution of the second term, followed by the integration over dipabd(pba results in the product of two integrals over 
two unconnected graphs obtained by eliminating the connecting edge {a, b}. 

Now we are in the position to describe a 
Proof of Theorem |II.7[ Substituting Eq. which expresses the higher order vertex functions in terms of the 



respective edge terms, into Eq. m, we obtain the following expression for the partition function in the form of a sum 
over the combined variables (7, ^ 

«(7:C) 1(C) 

^wBG(vr)=^(-i)^(^-«) n Kc,(7,o), Kc^)=n<'^u,- (9) 

(i,i) i=i 

Eq. is derived by means of re-grouping the factors correspondent to the lowest weights wic^ and keeping together 
terms correspond to the same immersed orbits Cj(7, ^) participating in the decomposition of 7 determined by ^. The 
sign factor in Eq. ^ is given by the product of the respective signs associated with the nodes of 7, and it is thus 



determined by the total number of intersections, N{'-f,^), introduced in Definition (II.5I. 
A similar to Eq. ^ expansion for the partition function of the Fermion model is 

Zg3(^,^;W^)- ^ n (-^(Q(7,0)) n (10) 

(7,«) J = l J = l 

It can be rationalized as follows. We represent the exponential under the average in Eq. Q as a product of the vertex 
exponentials (labeled by a) and expand the vertex exponentials in the natural bilinear combinations of the Grassmann 
variables ip. Each term of the expansion is naturally labeled by a set ^ of partition variables. We further compute the 
expectation value of each individual term in the expansion, using the Wick's theorem. Since the two-point correlation 
functions of the Grassmann variable is, {^Pab^ba) = (^ab, when {a, 6} € Qi and, {^Pab^ba) = 0, otherwise, the set ^ of 
the partition variables provides a non-zero contribution to the partition function if and only if it satisfies the following 
property: If {a, 6} participates in the local partition at node a, then {6, a} = {a, b} participates at the local partition 
at node b. For a partition ^ that satisfies this property we can build the associated cycle 7 consisting of edges {a, b} 
that participate in the local partitions. We further re-group the factors ipab'ij,c''Wbc'' ^"-c by keeping together the terms 
that correspond to the same immersed orbits Cj{'j, ^) participating in the decomposition of 7. The decomposition is 
determined by ^, followed by evaluating the expectation values using the Wick's theorem, in particular with making 
use of the form of the two-point correlation functions. After re-grouping the fact ors correspondent to each immersed 



orbit, and combining together the Wat\ Gab and <;ac'' factors, this results in Eq. (llOl) 



By Definition |II.6 the sign factors in all the individual contributions to the partition functions Zwbg(W^) and 



Zq3(^,(t;VF) (labeled by respective suitable orientation, (7,^)) arc the same, which proves the statement of the 
theorem. | 



IV. CONSTRUCTION OF A SUITABLE ORIENTATION 



In this Section we show how to construct efficiently a suitable orientation, required to compute the partition function 
of the WGB model efficiently, according to the Theorem II. 7 The construction of the suitable orientation is split 
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FIG. 2: Construction of the extended graph associated with Q. A degree three graph Q is constructed on an intermediate 
step. Q is obtained from Q by extending vertices. is obtained from Q by replacing linear segments with elementary edges. 



into two steps. First step, discussed in Section [TV A[ describes transformation of the original irreducible graph to the 
respective extended graph, Q'^, which is a homogeneous graph of degree three with even number of vertices. We also 
show in this Subsection how a Kasteleyn orientation of edges on the extended graph generates a so-called induced 
orientation on the original graph. Second step, undertaken in Section |IVB[ focuses on showing that an induced 
orientation of the original graph is always suitable, i.e. satisfies the condition of Theorem |II.7| 

Note that Proposition III.l reduces the problem of identifying suitable orientations to the case when Q is irreducible. 
Any irreducible graph of degree two is a cyclic graph, for which suitable orientations are identified in an obvious way. 
Therefore, throughout this section the graph C/, on which our models are defined is assumed irreducible and of degree 
three or higher (even when not explicitly stated), which does not lead to any loss of generality. 



A. Extended Graph and Induced Orientation 

Construction of the extended graph, Q*^, illustrated in Fig. ([2]), is done in two steps. We start from an original 
irreducible graph Q of degree three or higher, which, in particular does not have vertices of degree one, and extend 
its vertices of degree three or higher to circles (or polygons) so that the adjacent edges of the original graph are 
connected to the polygon vertices. The procedure results in an irreducible graph Q of degree three, which is however 
not necessarily homogeneous degree three, as it may have linear segments consisting of vertexes of degree two. We 
further replace the linear segments of Q with elementary edges to obtain an irreducible homogeneous graph Q'' of 
degree three. 

Definition IV. 1 (Extended Graph. Extending Polygons.) The constructed above (in the preceding paragraph) 
irreducible homogeneous graph Q'^ of degree three is called the extended graph associated with an irreducible graph Q . 
The polygons that replace the vertices of degree three and higher are called the extending polygons. 

We continue working the irreducible degree three homogeneous degree three graph Q'^ and observe that by construc- 
tion it has an even number of vertices. Thus, following Kasteleyn [TSl [H] (see also [57] for a comprehensive review) 
we construct the so-called Kasteleyn edge orientation on the graph. (See also Fig. [3^ for an illustration.) 

Definition IV. 2 (Kasteleyn edge orientation.) An edge orientation of Q'^ (which is a planar graph with even 
number of vertexes) is called Kasteleyn edge orientation of the graph if for any ( even or odd) face of the graph (possibly 
under exception of the outer face) the number of clockwise- oriented (negative) edges is odd. This is formalized in terms 
of the vector, cr'^ = (tr^jj = ±l|{a,5} G G"^), where crj^^ = if the arrow/orientation assigned to the undirected edge 
{a, 6} goes from a to b, and crjjj, = — 1 otherwise. 

Furthermore, G"^ is also equipped with a natural triplet orientation, built using the following rule: we set ^^"'''^ = 1, 
if going from 6 to c represents a left most turn at a, and <rf°'''^ — —1 otherwise. (See Fig. [sjD for an illustration.) <; is 
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(b) 



FIG. 3: Fig. (a) illustrates construction of the Kasteleyn orientation of edges - the three examples (with one, three and five 
clockwise edges per face respectively) are all legitimate Kasteleyn orientations of the selected internal face. Fig. (b) illustrates 
construction of a triplet orientation of a vertex of degree three (i.e. a vertex of Q'^). 




(a) 



(b) 



(c) 



FIG. 4: Fig. (a) shows two examples of possible choice of edge and triplet orientations within a linear segment of Q. The two 
examples are compatible, i.e. they both results according to Eq. (Ill in the same orientation of the respective single edge in C/*^ 
replacing the linear segment in Q. The number of arrows, marking both edge and triplet orientations within the linear segment 
of G, directed against the resulting arrow/orientation of the respective edge of Q"^ is even. Fig. (b) shows fragment of complete 
Kasteleyn orientation on and the related induced orientation of the respective triplet on G. Fig. (c) illustrates how two 
equivalent ways of counting orientations on Q"^ results in the same triplet orientation on Q, correspondent to the triplet from 
Fig. (b): the total number of the upper/left orientation within either of the paths (shown on left and right figure) which are 
against the resulting (red) arrow on Q is even. 



a vector constructed from all the triplet orientations of Q*^. The apparently sloppy construction is actually rigorous, 
since a planar graph, in particular Q'^, is supplied with a cyclic (e.g., counterclockwise) ordering of edges, adjacent to 
any node. 

Definition IV. 3 (Kasteleyn (Complete) Orientation.) A constructed above triplet orientation is called the 
left triplet orientation of . A (complete) orientation ((T^,<?^) of an irreducible homogeneous graph Q'^ of degree 
three, represented by a Kasteleyn edge orientation (J^ and the left triplet orientation ^ is called a Kasteleyn ( complete ) 
orientation of Q*^ . 

Inducing orientations. Our next task becomes to generate from (cr"^, ^'^), defined on Q'^\ a set of useful orientations 
(cr, (;) on the original graph Q. This task is solved in a number of steps. We first build a set of orientations (ct, ^) on 
Q by staying with the left triplet orientation for vertices of degree three and choosing the edge and triplet orientations 
on the linear segments of Q (end vertices excluded) in an arbitrary way with the only constraint per linear segment 
that the product, 

of all the edge an triplet sign factors reproduce the sign factor for the corresponding edge of G'^- See Fig. for 
illustration of the construction. Then, we turn to constructing orientation on Q, first building respective components 
of (<t,^) for degree two vertices to be identical to the respective (and just constructed) orientations of the degree two 
vertexes on Q. The components of <i related to degree three and higher vertices of Q are defined as the products of the 
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FIG. 5: Illustration of the lifting procedure. Top left: an example of Q (black) and an immersed orbit C (red). Three other 
sub-figures show three different (but equivalent in the sense of their respective ecr=,4=(C"^) contributions, see Proposition IV. 5 1 
lifted version of C on the respective extended graph, Q"^ . 



edge and inner triplet orientation factors of cr*^ respectively, over the paths that connect the corresponding edges 
going over the extending polygons in either (counterclockwise or clockwise) direction, as shown in Fig. [4jD,c. The 
Kasteleyn nature of the complete orientation (cr'^,^'^) ensures the independence of <; on the choice of the directions of 
the connecting paths. 

Definition IV. 4 (Induced Orientation.) A (complete) orientation (cr , <;) of an irreducible graph, constructed above 
(in the preceding paragraph) is called an induced (complete) orientation, or more specifically an orientation induced 
from a Kasteleyn edge orientation cr'^ of Q'^ . 

Lifting immersed orbits to the extended graph. Induced orientations are constructed in a way that the 
e(C)-factors for the immersed orbits in Q, defined by Eq. ([6]), can be studied by lifting them to the extended graph. 
Given an immersed orbit in Q we first lift it to Q by connecting the ends of the linear segments of C with paths going 
over the extending polygons. By further replacing the linear segments with elementary edges we obtain a lifted orbit 
in Q'^ . Obviously, by construction a lifted orbit, , is not unique and is specified by the choice of the directions 
(counterclockwise/clockwise) of the connecting paths. Lifting procedure is illustrated in Fig. [sj 

Proposition IV. 5 Let (cr, ^) on Q be induced from a Kasteleyn orientation (cr'^,^'^) on Q'^ , and be any lift of an 
immersed orbit C of Q to the extended graph Q'^ . Then Eo-s.^^ (C*^) — Sa-,^{C). 

Proof: The statement follows directly from the construction of the induced orientation. For £rT.q{C) the products 
of edge an inner triplet factors over the linear segments reproduce the edge factors cr'^ of the corresponding edges of 
Q''. The triplet factors of ^ for vertices of degree three or higher reproduce the products of the edge and inner triplet 
factors of (cr'^,^'^) over the connecting paths. The triplets that connect the linear segments to the connecting paths 
always come in pairs with two right/left turns for the counterclockwise/clockwise choice of the connecting path. | 

The important role of the Kasteleyn complete orientations and the corresponding induced complete orientations is 
described by the following statements. 

Proposition IV. 6 A Kasteleyn (complete) orientation of an irreducible homogeneous planar graph of degree three, 



Q'^ is suitable. According to Definition \II. 6 this is obviously equivalent to the statement that s{C'^) = —1 for any 
embedded curve (simple cycle) C ^. 



Proof: The simple cycle partitions the plane into two connected regions, interior and exterior respectively. There- 
fore, = dA4 is the boundary of the "inside" connected component (the one whose closure is compact), represented 
by a planar graph A4 d Q"^ (the boundary is included). It is enough to prove the statement for the counterclock- 
wise orientation of C"^, throughout this proof all involved orientations are counterclockwise. We further represent 
e(C"^) = eCT=£<;<= as a product of edge and triplet contributions along the embedded curve C"^, and denote by hq, ni, n2 
and I the number of vertices, edges, and faces of A4 and the length (number of edges) of C^, respectively. Considering 



^ We may consider any irreducible homogeneous graph of degree three here, and not necessarily the graph Q'^ extended from an arbitrary 
irreducible planar graph. However, since the proposition will be applied in the following solely to the extended graph, Q"^, we choose 
not to introduce new notations for an arbitrary irreducible homogeneous graph of degree three, using Q'^ for this purposes here. 
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FIG. 6: Illustration of the e^^ computation for an imbedded orbit of a homogeneous degree three graph, used in the proof of 
the Proposition |IV.6"1 Arrow along the marked (red) embedded orbit indicates anti-clockwise (positive) direction of the walk. 
± , standing next to vertexes along the embedded orbit (simple cycle), correspond to the values of respective <; contributions. 
The number of edges adjusted to the ± vertex and belonging to the interior of the cycle is two/three respectively. 



the product of Ea-^idM) = — 1 (the equality is due to the Kasteleyn nature of cr'^) over all faces of Ai we arrive 
at (—1)"'^ = (— 1)"^^'£ct= (C^), since each internal edge participates in the product exactly twice and with opposite 
orientations. Therefore, £cr=(C^) = (— l)"i~"2-i^ 'pj-^g triplet contribution is estimated e^e(C"^) = (— l)"'^ with nn 
being the number of right turns of C^. Since each internal node of M is of degree three, whereas a boundary node 
a G Co has degree two/three (in A^) when C turns left/right at a, and (^("^^^^ — ±1, respectively (see Fig. ^ for an 
illustration). In view of the above argument we have for the number of edges 2ni = 3(no — I) + 3n/j + 2(1 — nji). 
This results in e^e(C"=) = (-l)"""'. Comb ining the obtained results for the edge and triplet contribution we arrive at 
e(C"^) ~ (_i)"o-ni+n2 _ with the second equality due to the Euler characteristic relation, no — ni +n2 = 1. | 

Corollary IV. 7 If an immersed orbit C in Q is lifted to an embedded orbit in Q'^ , then Sa-,<;{C) — —1, provided {cr,(;) 
is induced. 

Proof: Let C"^ be any lift that satisfies the condition of the corollary. Then eo-.^lC*) — Ect^.^" (C"^) = —1 by Proposi- 
tions HVj and |IV31 



B. Induced Orientations Are Suitable 



The following statement explicitly provides a broad class of suitable orientations on an irreducible graph. 
Lemma IV. 8 Any induced (complete) orientation on an irreducible graph is suitable. 

Proof: We are not presenting a complete detailed proof of Lemma [IV.8[ since it contains a large number of straight- 
forward verification. We rather provide a sequence of steps that leads to a complete proof with brief explanations of 
these verifications involved on each step. The main idea is to demonstrate that for an induced orientation the binary 



function 9(7,^), defined in the Definition II.6 does not depend on the partition data ^ (Steps 1 and 2) and then 
explicitly compute it for some special choices of ^ (Step 3). 

Step 1. We start with noting that the triplet component <; of any induced orientation satisfies the following set of 
local constraints. Let a be an arbitrary vertex of degree 4 or higher, and (^o, 61, ^2, ^3) be a set of its four arbitrary 
distinct neighbors (i.e., bj ^ a), cyclically ordered in the counterclockwise direction with an arbitrary choice of the first 
neighbor. Then, the following relation, between the pair-wise components of the left triplet orientations associated 



with a and (601^17^2,^3) (see Definition IV.3), holds 

(a) (a) _ (a) (a) 

This can be verified directly. See Fig. [7] for an illustration. 

Step 2. We further introduce a set of local transformations of the decomposition data ^, illustrated in Fig. [Tj and 
show that 9(7;^) stays invariant under these transformations. Modulo 2 the change in N{'^,^) is totally determined 
by the change of the intersection index of the edges involved in an elementary transformation. The change in the 
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FIG. 7: Example of the elementary transformation which eliminates an intersection of the segments Ifbn .by} and I{bi,b3} 
shown in red. The products of the triplet orientation factors involved remain the same according to Eq. ( |12[ ), while the total 
intersection number of the involved segments with any other (e.g. these shown in gray) remains unchanged modulo 2. The 
change in the number of components, n{j,^), by one is compensated by the identical (modulo 2) change in the intersection 
index, iVa(7, 




FIG. 8: A single vertex element of an even generalized loop (Z2-cycle) 7 of 5 is shown in red at the top/left subfigure. 
Respected elements of the preferred decomposition data £, (containing no intersections) lifted to Q"^ is shown in red at the 
bottom/right subfigure. 



product of the edge and triplet factors involved in the product of all e-factors in Eq. ([s]) boils down to the change 
in the product of two triplet factors involved in the elementary transformation. The change in the number 71(7, ^) 
of the decomposition components can be also analyzed locally, as illustrated in Fig. [7) Therefore, the invariance of 



(_1)9(7,€) under an elementary transformation can be verified locally and boils down to the validity of Eq. (12). 

Finally, we note that by performing the described elementary transformation each decomposition data ^ can be 
transformed to a preferred one, shown in Fig. |8j where for each a E j the adjacent edges are paired with their 
nearest neighbors (counterclockwise or clockwise) on 7. Obviously, any preferred decomposition satisfies the property 
Nah, I) = for any a G 7. 

Step 3. The previous step has reduced the proof to the case when the partition data of a decomposed cycle (7, ^) is 
of a preferred form. We further lift the decomposition components Cj(7,^) to Q'^ by choosing the connection paths on 
the extending polygons in a way that they do not intersect with each other as shown in Fig.[8j which is always possible 
for our special choice of ^. Therefore the immersed orbits Cj(7,^) satisfy the conditions of Corollary IV. 7 which 



results in 9(7,^) = (— l)^^''"^-'. Also, as noted above, we have -/Va(7,€) = for all nodes a G 7, so that A^(7,^) = 0. 
I 



We would like to emphasize that Lemma |IV.8 not only proves the existence of suitable orientations, but rather 



provides a practical way of building them by reducing the problem to finding Kasteleyn edge orientations on the 
extended graph. 
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GAUGE TRANSFORMATIONS 



Gauge transformation approach discussed in [4l [5l [28] , keeping the partition function invariant and the graph intact 
but modifying the factor functions, also apphes to the WBG modeL On the other hand, and as shown above, the WBG 
model is det easy, and thus any binary model derived from WBG models via application of a gauge transformation 
is also easy. This explains the logic of the main part of this Section — to describe the extended family of binary 



models reducible via gauge transformations to WBG models. Section V A serves for a brief review of the gauge 
transformations. Sections |V B|V C|V D| describe the gauge transformations reducing the three noticeable examples of 
dimer, ice and Ising models to respective WBG models. 



Section V E describes a complementary example of an #X-matching det-tractable model, introduced in , which 
cannot be reduced (to the best of our knowledge) via a standard binary Gauge transformation to a WBG model, but 
it is still reducible to a WBG model on a contracted graph. The reduction to an easy model is shown here with the 
help of generalization to intermediate model defined over 3-ary (non-binary) alphabet. 

A. Gauge Transformations for general binary model. Reminder. 

Let us first define a general Edge-Binary Graphical (EBG) model on Q with the following partition function 

Zebg — 

E n (13) 

where tt = [iiab = T^ba = 0, l|{a, 6} e Gi), and Va e Go, tt^ = {iTablb G Go;{a,b} e Gi) and fai-^a) are the cost 
functions associated with vertices. 

As shown and discussed in [H [S] HH] the partition function of the general EBG model ([Is]) is invariant under the 
set of linear, so-called Gauge transformations. 



ao 



(14) 



where the newly introduced (two per edge of the graph) Gauge matrices satisfy the following skew orthogonality 
conditions 

V{a, b}eGi: E ^-''(^' ^')Gba{7T, tt") = ,5(7r', tt"). (15) 

TT 

With the gauge transformation applied, the partition function of the model becomes 

ZeBG=J2II fM=J2 n \T.[U^-b^''-b,7rUyj faiT^aU ■ (16) 



The expression in the middle of Eq. ( 16 1 may be considered as a new graphical model. Notice that if all gauge-matrixes 



are non-singular (which is the case discussed in this manuscript) transformation from ( 13 1 to ( 16 1 is invertible, i.e 



for any set of G-transformation one naturally defines the reverse G ^ , bringing us back to the original formulation of 



Eq. (16) 



In 4-6j 121] the Gauge transformation was discussed in the context of the so-called Belief Propagation (BP)- 
gauge and the resulting Loop Calculus/Series, where BP constraints were imposed additionally to the general (and 



always maintained) skew-orthogonality constraints (15). The main use of the gauge transformations in the preceding 
publications consisted in fixing the gauge freedom even further according to the so-called Belief Propagation (BP) 
conditions. The BP equations, enforcing the respective gauge fixing conditions correspond to fixed points of the Belief 
Propagation message-passing algorithm popular in physics, computer science and information theory. This special. 



BP gauge, reduces the number of terms on the rhs of Eq. ( 16 ) requiring that only terms correspondent to tt = 1-edges 



forming a generalized loop (i.e. subgraph with all vertexes of degree at least two). The resulting series for partition 
function consisting of generalized loop only, is called Loop Series. 

In the following we show how some set of binary models on the planar graph, which are known to be easy from 
previous studies, are reduced to WBG forms under respective gauge transformation. 
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B. Example: Dimer model as a WBG model 

Gauge transformation allows to restate the original model in a new basis, parameterized by G, which is generally 
a useful freedom to use. To illustrate its general utility, we discuss how applying a gauge transformation allows to 
restate the dimer model as a WBG model on the same graph. First, let us restate DM in the general EBG form 

f^'-\^a) = n * { I: ^oThe^^isi (1^) 

We will build here a gauge transformation around a valid dimer (perfect matching) configuration, p. Then introduce 
the following 7r-dependent gauge (2x2 matrix with elements GabiT^ab,T^'ab) parametrically dependent on p) 

Gab = 1 ) Pab = Pba = 1, Or clsC Gaf, = q ^ ^ ■ (18) 

Obviously this gauge satisfies the skew-orthogonality condition ( [l5| . Furthermore, one finds that 

r 1. Ec^a^ac^O 

/i (^aSPa) = < V'^o.c/Wab, Pab = T^ab = T^ac = 1, & Y^dta'^'^ad^O, (19) 

[ 0, otherwise. 

Restating it in words, for the new representation an allowed (nonzero) configuration at any vertex corresponds to zeros 
on adjusting edges (there are two of them), of which one is necessarily "active", i.e. present in the perfect matching 



and thus taking nab = 1 value. After collecting the Wab pre-factors on the rhs of Eq. ( 19 1 into one common multiplier, 
one observes that this modified model is in fact a particular Wick model ([2|. The factor functions of the model is 
identical zero for configurations with at least two adjusted edges taking "active" tt = 1 value. (I.e. in this particular 
case for any G in ([2| all W'^'''^ with Sb{C) > 2 are simply zero.) Notice also that the WBG representation for the dimer 
model is not unique as it can be build around any valid dimer configurations. Evidently any two valid representations, 
described by the gauges Gi and G2 can be transformed to each other via respective gauge transformations, G2G^^. 

We conclude this Subsection mentioning briefiy yet another alternative reduction of the dimer model to a WBG 
model, utilizing a topological transformation (specifically a contraction) of the original graph Q. Consider a valid 
dimer configuration po as a set of edges {a, b} € Gi and contract each edge {a, b} G po to a point. This results in a 
new graph Q, whose nodes Qo — Pq are represented by the edges that belong to the reference dimer configurations, 
and whose edges d = Qi \po are represented by the rest of the edges of Q. Each valid dimer configuration p on Q 
generates a configuration tt of a binary edge model on Q that has a form tTq, = 1 if a S p, and tTq = 0, otherwise. (Here 
we consider an edge a £ Gi ^ Gi\po Gi of G as an edge of G-) It is easy to see that a vertex function fa{T^a) for 
a € ^0 is nonzero only when all components of tt^ are zeros (the case a € po), or when exactly two components have 
the value one (otherwise), and also that the local weights satisfy the conditions of the Wick binary model expressed 
in Eq. ^. 



C. Example: Ice Model for graphs of degree three 

The ice model is defined in terms of orientations of edges on the graph. A valid configuration/orientation is such 
that no vertex has adjusted edges all oriented onwards or inwards. Here we consider the case when all vertexes on 
G are of degree three. (Without loss of generality our consideration here can be extended to vertexes of degree not 
larger than three.) Weight of any allowed orientation of the graph is unity. This model, #PL — 3 — NAE — ICE 
according to the complexity theory classification, was discussed in [U [3] in lieu of its reduction to dimer model via a 
holographic algorithm. Related models were studied extensively in the mathematical physics literature P^24j . 

The ice model can be conveniently restated in terms of normal binary variables ^. This will require a graphical 
transformation consisting in breaking every edge, {a, 6}, by inserting a new auxiliary vertex a — b. Then, binary 
variables assigned to the pair of new edges are T:a,a~b — 0, nt^a-b = 1 if the direction of the arrow (for the original 
orientation of the edge) was a ^ b and T^aja-b = Ij '^b,a-b = for the b ^ a arrow/orientation. The partition function 



An alternative derivation, not requiring a graphical transformation, is discussed in Appendix 
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of the respective binary model becomes 



b,Trb,a-b)\, (20) 



r f / \ J Ij 3 6, C ^ a, S.t. TTa.a—b 7^ '"'0,0— c /rji \ 

Ja[T^a) - \ 0, Otherwise ' ^ > 



ga-^b{T^a) ~ I 



1 T^a,a-b 7^ T^b,a~b 

otherwise 



(22) 



where tt' = tt U (7ra,a-6 = 0, l|{a, 6} e t/i) and tt^ = (7ra_a-b = 0, l\b (5e(a)). 

Let us introduce the foUowing gauge transformation on aU edges of the newly formed extended graph 

G^:ib - ( -1 ' ^''^ 

Factor functions of the ice model in the new basis becomes 



3 



/a(7ra,a-l, 7rQ^Q-2, TTa^a-a) = ^ * N -1/3, Z^i ""a.a-i = 2 , (24) 



%/2 



0, otherwise 



(ja^b^T^'a) = { '^a,a-b = '^b,a-b = 1 , (25) 

0, otherwise 

where (a, a — i) with j = 1, 2, 3 marks three edges adjusted to vertex a. Accumulating the product of the overall 3/\/2 
terms into a common (normalization) multiplier, one immediately finds that the resulting model is of the Wick's type 
([2]). Indeed, given that Va G 5 5a(^^) < 4, one simply does not have any M^f") terms in the ice model version of ^ 
with Sa{G) > 4, and the only nontrivial vertexes are these with Sa{G) = 2. Notice, that this Wick feature would 
not apply to the ice model on planar graphs with vertexes of degree higher than three — one can still use the gauge 
transformation but the resulting transformed model will not maintain Wick's property. 

Note that transformation of the ice model, identical to the one discussed above (e.g. utilizing some alternative 
terminology) was also used as a show case of the holographic algorithms of Valiant in [IH5] . 



D. Example: Ising Model 

We consider the classical Ising model without magnetic field. Original formulation is in terms of binary variables 
associated with vertexes, however the following EBG representation is more appropriate for our purposes 



^/=X!( II I II 9ab{T^a,a-b,Trb,a 

n' KaeGo I \{a.b}eQi 



-b) I , (26) 

r I \ / Ij V&, ce(5g(a) T^a.a-b — T^a^a-c^ /f,-\ 

Uy^) - I 0, Otherwise ' 

n ( ^ ^ J ^' '^a,a—b '^b,a—b foa\ 

9ab['^a,a-b,'^b,a-b) ^ \ / , (2») 

'^a,a—b 7- '^b,a—b 

where we adopted notations of the previous Subsection for new (auxiliary) vertices breaking each edge of the original 
edge in two. 



Applying the Gauge transformation ( 23 ) , introduced above for the Ice model, to the Ising model one arrives at the 



following factor functions in the new basis 



1, 7ra,a-6 — nb,a-b = 

Jab{Tra,a^b,Trb,a~b) = h + ^J■) * { {l - fJ-) / il + ^J-) , '^a,a~b = ■^b,a-b ^ ^ ■ (30) 

0, otherwise 



k = 2 



2x 
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2 + l = 3 = 3!!^-(2jt-l)!: 



yt = 3 




3 + 6 + 3 + 2 + l = 15 = 5!!^(2/t-l)!! 



FIG. 9: Illustration clarifying the Wick-feature of Eq. (291 



With the product of all the pre-factor terms, 2~l''6'(''^l/^"'"^ and (7 + accounted for in the overall multiplier, 
the resulting model is one of the WBG models of ([2|. This Wick feature of the modified graphical model is ob- 
vious for degree two vertexes of Eq. ( 30 ) , and the only part which requires additional clarification concerns the 
relation between a configuration of Eq. ( 29 ) with X]&~a "^afi-a > 2 and its pairwise decomposition. Indeed, for a 
^2-cycle 7 and a vertex a € 70 of valence 6-y{a) = 2k with respect to 7, there are (2k — 1)!! local partitions with 
{{2k — 1)!! + 1) /2 and {(2k — 1)!! — 1) /2 of them having modulo two intersection index equal to zero and one, respec- 
tively. Therefore according to the Wick rules, the cumulative local contribution associated with the = 2k vertex is, 
{{2k - 1)!! + l)/2- ((2k - 1)!! - 1) /2 = 1, that is fully consistent with Eq. ([29]). This simple combinatorial relation 
is illustrated in Fig. |9] for k — 2 and fc = 3. 



E. #X-matching model 

The #X-matching model is defined as follows [1H3] : On a planar bipartite graph, where Vi and V2 are bi-partitions 
of the full set of nodes and the nodes in Vi all have degree 2, consider matchings (monomer-dimer configurations) 
weighted in a way that all dimers (edges) have arbitrary weights, monomers on Vi are all of unit weight and the 
weights of the monomers on V2 are equal to minus the sum of the dimer weights of edges emanated from the vertex. 

It is convenient to view the bipartite graph in the original definition of the #X-matching model as a graph with 
the the sets of nodes V and edges E represented hy V = V2 and E — Vi, respectively. We will adopt notations 
similar to these used in Section V C[ devoted to the Ice model, calling a, 6, • • • the vertexes of 1^ = V2 and a — b, 



vertexes belonging to E = Vi. Then the partition function Z^x of the ^^X-niatching model has the following EBG 
representation on the graph Q — (Go,Gi) — (V,E) 



(31) 



%X I II /a(^") I n 9ab{Tra,a-b,T^b,a-b) , 

J2br^a^ab'^CL,a-b, J2br^a'^a,a-b ^ I , (32) 

0, otherwise 

9ab(^a,a-b,l^b..a-b) - | ^ ^ g ^ q , (33) 

Our task in this Subsection is to show via a gauge-transformation approach, alternative to the match-gate approach 
of [IHS], that the partition function of the #X-matching model, defined as the weighted number of the monomer- 
dimer configurations, is det-easy. To achieve this goal we will need to introduce, as in the above Subsection, a gauge 
transformation reducing the model to another model known to be det-easy. However, the gauge transformation of 
this Subsection is principally different from these discussed in the rest of the paper, as here we will need to consider 
generalization from binary-setting to a g-ary, with g = 3, alphabet. For this purpose we restate Eqs. ( 3l][33 ) as the 
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following 3-ary Edge Graphical Model over the set of variables tt = {iTa.a-b} that attain the values TTa.a-b = 0,±1: 

Z^X ( n /'^(^a) ) n 9ab{TTa.a^b,T^bM~b) , (34) 

{1, ybr^ a TTa^a^b = 0, 

WabTra,a^b, TTa.a-b = ±1 and Vc ^ 6 TTa^a-c = , (35) 
0, otherwise 

— /- ~ \ JO, '^a,a — b — ^b,a—b — 1 /'o^?^ 

9ab{7ra^a^b,7r,^a^b) ^ < , otherwise • 



To see equivalence between the binary and q-ary representations let us walk backwards, from Eqs. ( 3l][33 1 to Eqs. (34 



hen 



36), and associate with any configuration tt of our 3-ary model a binary configuration tt by setting Ha^a-b = 1 wf 
T^a,a-b — 1, and Ha.a-b — whcn T^a.a-b = 0, —1. Then, we do partial summation over all ifa^a-b = 0, —1 correspondent 
to TTa^a-b = 0. In a general case the result of such partial summation would be non-local with respect to the graph. 
However, in our special case the resulting model happens to be graph-local. Indeed, the weights gab(J^a,a-bTT^b.a-b) in 



Eq. (36 1 actually depend on the reduced variables na,a-b and nt,^a-b and the partial summation can be performed for 
each node a £ Qo independently. Finally, we arrive at the binary formulation of Eqs. ( 3l]|33 ). 



One may wonder why changing from binary to 3-ary representation is useful? In fact it is very useful because of 
a greater freedom of further transformation the 3-ary representation offers. Indeed, we will see below that relating 
TT variables to tt variables in a new way allows to reduce the 3-ary model to the already known det-easy model 



(specifically, a dimer model, already shown reducible to a WBG model in Section VB). 

Namely, we set iTa^a-b = 1 when TTa,a-b — il, and iTa^a-b = when i^a.a-b = 0. The resulting binary model turns 
out to be the dimer model on the reduced version of Q (no half-edges only full edges) with the dimer weights —WabWha 
for the {a, 6} edges. This follows from the following arguments. Consider a link {a,b}. If 'Ka,a-b = 1 and ■Kb,a~-b = 
we have TTa,a-c — for all c ~ a with c 7^ 6, and also gab{T^a,a-bTT^b,a-b) — 1, so that the partial summation over 
T^a.a—b — il Can be performed independently, which results in 0. If TTa,a—b — '^b,a—b = 1, we have TTa.a—c 

= for 

all c ^ a with c ^ h, and also T^b^b-c — for all c ~ 6 with c ^ a. This allows the summations over iTa^a-b = il 
and iTb,a-b = ±1 to be performed independently, thus resulting in the factors of for T:a,a-b = '^b,a-b = 1, and 
WabWbaT^a,a-bT^b,a-b, Otherwise. QED. 

It is worth to note that in terms of the language of loop towers for g-ary models, developed in our earlier work on 
loop calculus for g-ary alphabets |28j . the partial re-summation based on partitioning the alphabet into one special 
letter (tt = 0) and the rest (tt = ±1) should be viewed as a partial gauge fixing that corresponds to the lowest level 
of the loop tower. 

VI. DISCUSSION AND FUTURE PLANS 

In the present manuscript we have identified a broad class of models that allow for an easy solution on planar 
graphs. Our approach is based on the gauge invariance of vertex models on graphs that has been implemented in our 
previous work in the context of Belief Propagation and Loop Series and the notion of a Wick vertex model, referred 
to as WBG. The WBG models are described in terms of an arbitrary choice of weights associated with edges of the 
graph, and explicit construction of the high-degree (even coloring) vertex weights from the respective edge weights 
according to the so-called Wick rules. We have established an equivalence between the WBG models and th e Ga ussian 



fermion (Grassmann variables based) models on the same graph, referred to as the VG'^models (Theorem II.7I. The 
models are "easy" as their partition functions are equal to Pfaffians of square matrixes of the size equal to twice the 
number of edges of the graph. We have also demonstrated that the dimer model that constitutes the "solving tool" 
of the holographic algorithm of Valiant [iHS] is gauge equivalent to a WBG model. The other two components of the 
holographic approach include graphical transformations (extensions of the original graph) and linear transformations, 
which can be described as particular cases of gauge transformations using the language adopted in this manuscript. If 
one represents the partition function of the dimer model on the extended graph, obtained by using the gadgets of the 
holographic algorithm, as a Gaussian Grassmann integral and integrate over the variables related to the new vertices 
that result from the graph extension, the result will be represented by a Gaussian Grassman integral on the original 
graph. Therefore, our approach provides an explicit and invariant description of the models that can be treated using 
the holographic algorithm as the models gauge equivalent to WBG models. It also provides with an alternative way 
to obtain an easy solution, while sticking to the original graph all the time. 

Some problems we plan to address in the future, extending on the results/approach described in the manuscript, 
are: 
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Generalize the approach to account for det-easy planar graphical models defined over q-ary alphabet with q > 2. 



As illustrated in Subsection V E there exist an interesting new class of more general g-ary models which are 



reducible to det-easy binary graphical models. 

• Use the described hierarchy of easy planar models as a basis for efficient variational approximation of planar but 
difficult problems, and possibly more generally non-planar models. Then, on the next level, build a controlled 
perturbative corrections to the variational result. Notice, that this approach may also be useful for building 
efficient variational matrix-product state wave functions for quantum planar models, e.g. in the spirit of [22]. 

• Analysis of Wick Gaussian models on surface graphs of nonzero genus. This work will extend the classical 
results [m [3OH35I stating that partition function of dimcr models on graphs embedded in surface of genus g, 
so-called surface graphs, is expressed as a sum over 2^^ Pfaffians, each correspondent to a "spinor" structure 
parameterizing the equivalence classes of Kasteleyn orientations on the locally planar graphs. The approach 
developed in this manuscript allows for a relatively straightforward extension to the case of surface graphs. In 
particular, we plan to show in a forthcoming publication that Kasteleyn orientations on the extended graph 
generate all 2^^ spinor structures on the embedding Riemann surface. 

• Study Wick Gaussian models on non-planar but Pfaffian orientable or fc-Pfaffian orientable graphs (thus any 
dimer model on surface graph of genus g is 2^f -Pfaffian orientable), also utilizing the new results from graph 
theory on the Pfaffian orientability • 

• For the case of a general EBG model we will apply the statistical variational principle to build the "best" 
approximation in a form of a WBG model and address the problem of finding efficient ways to account for the 
corrections. 
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Appendix A: Topology behind the door: Self-intersection invariant of immersions and spinor structures 

The main result of this manuscript is represented by Theorem [n.7[ The proof of the Theorem rests on a com- 
putation of the relevant Gaussian Grassmann integral using a certain expansion, combined with the Wick Theorem 



and Lemma IV. 8 In this Appendix we present a brief quahtative discussion (without proofs) of the topology that 
stands behind the presented combinatorial proofs presented in the main body of the manuscript. This will connect 
our results to the approach developed in [Ml [3S] for the dimer model on surface graphs, and also provide some ideas 
on how our results on the equivalence between the Wick and binary models can be extended to the surface graph case 
(the details of this description and related proofs will be published elsewhere) . 

A topological nature of the equivalence between the WBG and G'^ models can be seen by comparing Eqs. ^ 



and (10) that play a major role in the proof of Theorem II. 7 Each contribution in both sums can be represented 
as a product of the individual contributions labeled by the immersed orbits C'j{j,^). The individual contribution 
associated with an immersed orbit C has the same absolute value for both expansions, they differ by the sign factors 
only. For the G'^ model the sign associated with an immersed orbit C is — e(C) = (— l)«(C)+i (where the last expression 
defines s(C) as a function of e{C) introduced in the main text), in the WBG model case the sign is (— 1)^('^) with 
N{C) being the number of self-intersections of C. Note that the intersections between different orbits in the expansion 
of 7, although apparently affecting the sign factor in Eq. (|9|, can be actually eliminated from the consideration, since 
in the planar case two orbits always intersect at an even number of points. In this appendix we will argue that s(C) 
can be associated with the unique spinor structure on , and the equivalence between the models originate from the 
topological formula, N{C) + 1 = s{C) (mod 2), that relates the self-intersection number to the spinor structure. 

We proceed with more precise formulations. Note that the statement of Lemma [IV. 8[ can be interpreted in the 
following way: The function (7(7, ^) of 7 and ^ defined by Eq. ([8| can be viewed as a Z2-valued function defined 
on the sets of immersed orbits, represented by Cj{j,^), and it satisfies the property 9(7, ^) = 0. 9(7, ^) consists of 
three contributions: the number A^(7, ^) of intersections and self-intersections of the involved orbits (naturally modulo 
two), the number ^(7, ^) of the orbits involved, and the contributions s {Cj (7, ^)) that arise from the individual orbits, 
where s = (log(e)/(i7r). What wc intend to illustrate is that by combining the modulo two numbers of intersections 
and the modulo two number of orbits in a decomposition of a planar Z2-cycle 7 we arrive at a topological invariant 
p(7j — ^ili ^)+"'(7: €) G ^2, referred to as the self-intersection invariant, whereas the individual orbit contributions 
s {Cj{j,^)) and the corresponding sign factors e — {—ly are associated with a spinor structure on the plane M^, the 
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(a)An example of a smooth closed 
trajectory. Parametrization 
chosen for plotting is r = 
(nW,r2(i)) = (sin(t),sin(2t)), 
where t changes from to 27r. 



(b)Dircction of velocity (angle), 
ip{t) = arctan(2 cos(2t)/ cos(t)), for a particle 
moving along the immersed orbit from (A), 
thus r are coordinates of the particle in M^, as 

a function of "time" t. The invariant, 
measuring total number of angular rotations, is 
for this example. 



(c)Trajectory of the 
particle in the 
phase-space, x = 

(nW,r-2W,¥'W) 6 
M3 = R2 X 



FIG. 10: Illustration of an immersion in of an orbit with one crossing and its phase-space (particle) counterpart. 

original planar graph and its extension are imbedded into. Combining the self-intersection invariant with the spinor 
structure related contributions we arrive at the Z2-invariant q — p + s that depend on the planar Z2-cycle 7 only, and 
q = 0, since all cycles in are contractable. This will be done by interpreting the orbits as closed trajectories of a 
free particle in the plane, also referred to as a free planar particle. 

1. Immersions, self-intersection invariant and spinor structures on 



This Subsection introduces some useful objects that are new to the manuscript, as defined for continuous spaces, e.g. 
rather than for the graphical structure as in the main part of the manuscript. This continuous spaces approach 



will be used in the following Subsection for a topological illustration of Lemma IV. 8 The arguments presented 
there can be easily converted into a proof, however, this task goes beyond the scope of this manuscript. We will 
start with introducing a special type of smooth maps (trajectories) C : 5^ — M^, referred to as immersions and 
denoted C : 5^ M^. The main reason for bringing the immersions to our discussion is that the immersed closed 
walks, t hat re present the immersed orbits Cj(7,^) in the decomposition of a Z2-cycle 7 on a planar graph ^ C in 
can be deformed in a natural way to immersions Cj : S'^ S-> and studied using topological methods. 



Lemma 



IV. 8 



This rationalizes the terms immersed closed walk and immersed orbit, we have introduced in the context of graphs. We 
further discuss the self-intersection invariant p : Imm(5'^,M^) — >■ Z2 that associates with an immersion C : ^ M.'^ 
the number of self-intersections modulo two. The self-intersection invariant provides a topological interpretation of 
the intersection number iV(7, ^) that appears in Lemma IV. 8 In the context of the self-intersection invariant we will 
refer to Smale-Hirsch-Gromov (SHG) theorem |37| that reduces the problem of studying the topological properties of 
immersions C : ^ to a. much simpler problem of studying their phase-space counterparts /(C) : M.'^ x . 

In particular, the SHG theorem rationalizes a natural appearance of the phase space = M.'^ x . We proceed with 
introducing the notion of the spinor structure on viewed as a binary function s : LM^ — Z2 that associates with any 
phase-space trajectory C : S*^ — ?> a number s{C) e Z2 and can be also defined for the immersions C : 5^ by 

s(C) = s(/(C)). The spinor structure provides a topological interpretation of the factors e (^^ (7,^)) = (^—lyi^^^^'^'f) 
in terms of the value of the spinor structure on the relevant immersions. Finally, we briefly discuss a topological 
relation between the spinor structure and self- intersection invariant that stands behind the statement of Lemma |lV.8| 
We start with introducing a concept of an immersion. To that end we consider a particle moving smoothly 
over a closed trajectory in a plane two-dimensional space K^, where thus smooth planar trajectory C : 5^ — >■ 
parameterized by r{t) = {ri{t),r2{t) where the "time" t belongs to S^. The trajectory is called an immersion, which 
is denoted by C : S'^ if the velocity stays non-zero, i.e., r{t) ^ 0, at all t. In the following the key space will 

be the iso-energetic shell = x 5^ of a free planar particle so that for x — {ri{t), r2{t), 0{t)) e x , where 
the angular variable 6 G describes the velocity direction (i.e., the normalized particle velocity) r{t) 1^(^)1^ G S^. 
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Following the terminology, commonly accepted in mathematical physics, we will also refer to the iso-energetic shell as 
the particle phase space. For a smooth planar trajectory C with a non-zero at all times velocity, i.e., an immersion, 
we denote by /(C) its phase-space counterpart with the normalized velocity, which clearly forms a trajectory in 



(see Fig. 10 for an illustrative example). 

Our next step, as outlined above, is to describe the self-intersection invariant, defined as the modulo two number 
of self-intersections i{C) — N{C) (mod 2) of a smooth planar trajectory C : S"^ — > M^. It is intuitively clear that 
i{C) is a topological invariant only for immersions S-> M^, i.e., it stays the same only if we deform a trajectory 
so that it stays an immersion all the time, since in this case the self-intersections can appear and disappear in pairs 
only ^. Immersions can be conveniently studied using the SHG theorem that can be applied for our purposes in 
the following way. Interpret the procedure of associating with an immersion C its phase space counterpart f{C), as 
described above, as a map / : Imm(5'^,M^) — >■ LM^ from the space of immersion of a closed trajectory on a plane 
with everywhere nonzero velocity to the space LM^ of closed trajectories in — M."^ x S^. The SHG theorem, 
applied to our case, claims that / is a weak homotopy equivalence. In particular it implies that given two immersions, 
C, C" e Imm(5'^,K^), to answer the question whether C can be deformed to C", while staying an immersion during 
the whole deformation, is equivalent to finding out whether /(C) can be deformed to f{C') without any further 
restrictions. Note that the SHG theorem explicitly classifies the immersions ^ M^: since the coordinate space 
is contractible, only the velocity component 9c '■ ^ of the phase-space trajectory /(C) : S*^ ^ x S*^ 
matters, with 9c being fully homotopically determined by its degree deg{9c) G Z that describes the total number of 
rotations the velocity direction makes over the whole trajectory. In particular two immersions C, C £ Imm(S'^,M^) 
are equivalent if and only if deg{9c) = <i-eg{9c'), i.e. the direction of velocity (the angle) makes the same number of 
rotations over both trajectories. 

Consider concatenation of two trajectories Ci*C2 that share the starting point point in M^, understood as trajectory 
of a particle cycling C2 first and then going over the Ci. The concatenation results in 

i(Ci * C2) = i{Ci) + i{C2) + Ci-C2 + l, (Al) 

where Ci • C2 denotes the modulo two intersection index (the modulo two number of intersection points) of Ci and 



C2. Eq. (Al) reflects the fact that the number of self-intersections of a concatenation of Ci with C2 is the number 
of self-intersections of Ci plus the number of self-intersections of C2 plus the number of intersections of Ci with 
C2 plus one (the self-intersection at the point of concatenation), which, as opposed to the self- intersection index, is 
well-defined for any continuous closed trajectories, not necessarily immersions. Note that for the planar case discussed 
here Ci • C2 = (intersections appear and disappear in pairs), stated differently it reflects the fact that all cycles in 
are contractable. Nevertheless, we keep this contributions, since for our application it will be convenient to treat 
intersections and self-intersections on equal footing. By defining p{C) = i{C) -I- 1 we arrive at 

p(Ci*C2) =p(Ci)+p(C2) + Ci •C2. (A2) 

By the SHG theorem the topological classes of immersions ^ are fully described by the topological classes 
of their phase space counterparts /(C) : — )■ = x 5^, the latter being labeled by the first homology group 
Hi{M^) = Z, counting the number of rotations of the velocity orientation along the immersion. The self-intersection 
invariant, therefore, is also defined in terms of the natural factorization p : Z — ^ Z2 (i.e., it can be viewed as a 
composition LM^ Hi{M^) Z2 where the left map associates with a closed phase space trajectory the number 
of the velocity rotations), and satisfies the property 

pill +72) =p{li) +P{l2) +5(71) -5(72) =p(7i) +^(72), (A3) 

where (7(7) is a closed trajectory in that is obtained by reduction (projection) from by simply ignoring 
information about orientation of the particle velocity. (7(71) • 5(72) stands for the intersection index (mod 2) of the two 
reduced closed trajectories, and for the planar case considered in the paper it is equal to zero, i.e., (7(71) • 5(72) = 0. 
Consider a composite trajectory defined as a union of smooth closed and not necessarily intersecting trajectories 

(immersions). In this case we define p ^X]j=i Q;) = * (Llj=i ^i) + which generalizes the definition of a single 



If the condition of having everywhere non-zero velocity is relaxed, this would be not true anymore. This can be illustrated using a 
simple example. One can imagine a closed trajectory of the shape of "figure eight" from Fig. |10| being deformed in a way that one 
half of it becomes smaller turning into a point with a cusp, followed by smothering the cusp, which results into a circle. The number 
of self-intersections changes by one. Note that at the stage of the deformation when we have a cusp, the velocity of the trajectory at a 
cusp turns to zero. This means that the map C ; 5^ — > M-^ at this stage is not an immersion, which actually allows the parity of the 
self-intersection number to be changed. 
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FIG. 11: Schematic illustrations for details of the transformation procedure from an element of a directed orbit on a graph to 
its particle-trajectory smooth analog embedded into and consequently into . Notice that if the direction of the orbit is 
reversed the trajectory f{C') in will change. Underlying spinor structure is not shown on the figure. 



smooth trajectory (an immersion). Note that this definition ensures the property given by Eq. (A3 1. Summarizing, 
we have an invariant ^(7) defined on Hi{M^) = Z that satisfies Eq. (A3 1 and if a complex trajectory is decomposed 
into a set of smooth closed trajectories we derive 



Ep(^J-) +Y.9iCi) ■ 9{Ck) - 



(A4) 



i=i 



]<k 



This completes our discussion of the self-intersection invariant. 

We are now in a position to introduce the notion of a spinor structure, as outlined in the beginning if this subsection. 
A spinor structure can be interpreted as a map s : LAfi — > Z2 that associates with any phase-space trajectory 
(7 : S"^ — )■ X S'^ its parity s{C) = 0, 1. The spinor structure map s should satisfy the following two conditions: 
(a) Any continuous deformation of C does not change s(C), or more accurately (and generally) if C U C form a 
boundary of a two-dimensional domain then s(C) = s(C'); and (b) s{Cr) = 1 for Cr{t) — {r,t) being a phase space 
trajectory with the constant position r and the velocity making a full rotation. Obviously, the spinor structure can 
be also viewed as a map s : Imm(5^,M^) if we define s{C) — s(/(C)) for C : iS^ S-> M^. For this interpretation the 
condition (b) means s{Cr) = 1 where the immersion Cr represents a circle on K^, centered at r with a very small 
radius. Obviously, according to the Stokes theorem applied to the phase space, any spinor structure s(C) can be 
represented as 



siC) = s(/(C)) 



Aj{x)dx^ , 



f{C) 



(A5) 



where A is an Abelian curvature-free vector potential A, in the phase space — M."^ x . The curvature-free 



condition reads K 



0, thus expressing the condition (a) formulated above. ^From the definitions 



of the spinor structure and the vector potential, one finds that the value s(C) of the spinor structure map on any 
immersion C : S*^ is given by the modulo two total number of the particle velocity rotations along the 



closed trajectory (see Fig. 10: for an illustration). Note that for our planar case the conditions (a) and (b) define a 
unique function s, which means that there is a unique spinor structure on (this is not true in the surface case). 
Also, this unique spinor structure can be described by difi^erent gauge fields A that differ by gauge transformations 



djcj). This completes our discussion of the spinor structures. 



2. Topological interpretation of Lemma |IV.8l 



To demonstrate the equivalence of the combinatorial and topological definitions of e(C) we view a directed immersed 
orbit C" on Q*^, associated with an immersed orbit on as a continuous piece- wise smooth trajectory in M^. See Fig. 11 
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for a schematic illustration of the orbit-to-trajectory transformation. In very small neighborhoods of the nodes we 
deform C in an obvious way, so that we do not gain self-intersections in these small neighborhoods, to obtain an 



immersion C : S 



In the small neighborhoods of the nodes, where the deformations have been performed. 



the coordinate r, i.e., the component of the corresponding phase-space trajectory f{C') : S-^ — >■ is almost 
constant, whereas the velocity orientation 9 makes a finite rotation in the counterclockwise or clockwise directions for 
the left and right turns, respectively. Following f{C') on the deformation of C' back to C we obtain a continuous 
piece- wise smooth trajectory in that represents a orbit in M^, denoted with some minor abuse of notation 
f{C') : — > A'P, whose smooth pieces Bab and correspond to the edges and the nodes (velocity turns at the 
nodes) of C , respectively. The value of s {f{C')) can be consi dere d as a sum of the contributions from the smooth 
pieces, labeled by the edges and nodes of C", according to Eq. (A5). Indeed, we denote by Bab a smooth path in 



that represents the edge of in the embedding C M"^, whereas for a triplet a — ^ & — > c of we denote by B^'a} the 
shortest path on {rf,} x that connects the directions of the velocities at rf, on the paths Bab and Bbc^ respectively. 
Then the aforementioned edge and vertex contributions have the form 



Sa& = s{f{Bab)) 



Aj{x)dx^ , 



(A6) 



respectively. Naturally, e{f{C')) = [—ly^f'-^ )) is given by the product of the edge and node contributions and, 
provided the latter are represented by the Kasteleyn edge orientations and left triplet orientations, which can be 
always achieved by a gauge transformation of the field A, we reproduce the original combinatorial form, i.e., ec{C) = 
erifiC)), where the subscripts C and T stand for combinatorial and topological, respectively. 

To evaluate q{j) for a Z2 cycle we decompose it into n^j,^) immersed orbits Cj = Cj(7,^) and consider Cj and 



C" as continuous piece- wise smooth orbits in K^. We denote by Cj : S ^ 



e and C'J 



51 ^ 



the immersions. 



obtained by slight deformations of Cj and C" in small neighborhoods of the nodes of G and C/'^ , respectively, using 
the procedure described above in the context of C. It is easy to see that Cj, Cj, and Cj can be deformed to each 
other within the immersion space, in particular f{Cj), f{C'j), and f{Cj) represent homologically equivalent cycles in 
M^. Therefore, we have 



9(7) = 9 1 U = E + E ■ = E ^ (/(^^)) + E^^ (/(c-;)) +J2^3- c; 

U = l / j = l j<k j = l j = l j<k 



E ^ (m)) + - + E * (/(^;)) + E ^' • ^0 = E - (m)) +-+^^ 

J = l \J = 1 j<k I 3 = 1 



(A7) 



where N = N{"/,^) is the total number of intersections and self-intersections. Combining Eq. (A7) with th e rela tion 
^c{Cj) = Et {f{C'j)) = (— we complete our topological interpretation of the statement of Lemma IV. 8 



3. Topology on Surface Graphs: related approaches and future work 



A topological interpretation of Lemma IV. 8 discussed in this Appendix is based on establishing an equivalence 
between the definitions of e(C) in terms of a Kasteleyn orientation on the extended graph Eq. (|6| and a spinor 
structure on K^, respectively. This description allows generalization from to surfaces with finite genus. A 
relation between Kasteleyn orientations on surface graphs Q C and spinor structures on the surface has been 
established by Cimasoni and Reshetikhin [34] , who demonstrated that a combination of a Kasteleyn orientation and a 
valid dimer configuration on Q produces a spinor structure on M^. We have implemented a very similar, still different 
approach by associating a spinor structure on with a Kasteleyn orientation on the extended graph ^'^ C M^. In a 
forthcoming publication we will demonstrate that generalization of our construction solves the problem in the general 
surface case, thus allowing to relate partition function of a general discreet-variables Wick model defined on a 
graph embedded in a surface of genus g to a sum over 2^f contributions, each associated with a partition function of 
a fermion model of an allowed equivalence classes of spinor structures of Af ^ . 

It is also worth pointing out that Cimasoni and Reshetikhin |34| used the aforementioned self-intersection invariant 
N{C) implicitly since the quadratic form q they have used to study the dimer model on surface graphs and the 
one considered here (in our planar case q = 0) are actually the same. However, the objects that naturally arise 
in the studies of the dimer model on surface graphs are represented by simple unions of non-self-intersecting loops 
[3TJ [351 [331 [3S] , and, therefore, the issue of intersections is never raised in the context of the dimer model. On the 



25 



contrary, intersections play a key role in establishing the correspondence between the Wick and fermion models, and 
should be handled explicitly. To the best of our knowledge, in the context of 2d statistical mechanics for the first time 
the self-intersection invariant, of the type discussed above, however in a gauge of another form, has been brought up 
by Kac and Ward [lOj to come up with an exact solution for the two-dimensional Ising model on a square grid. In the 
later work, where the Majorana spinors on irregular lattices have been studied j3Hl [35] , the authors have been using 
the complex valued phase factors e** to establish relations between the free fermion and binary (Ising) models. This 
can be viewed as an extension of the approach of [10' from a regular to irregular lattices. In this manuscript we have, 
in a sense, followed j34. by working with the "Kasteleyn" phase factors ±1. The equivalence between the approaches 
that use different phase factors can be established using a gauge transformation of the field A that describes a spinor 
structure according to Eq. ( A5 ) . 



Appendix B: Extending the gauge group 



In this Appendix we discuss an alternative way of reducing the three illustrative models, dimer, ice, and Ising 
discussed in Section |v] to the WBG model described by Eqs. (2j3). In Subsection VA the dimer model has been 
converted to the WBG model via an example of the gauge transformation described in Eqs. ( 14|15|16" ). The ice and 
Ising models strictly speaking do not belong to the class of vertex models, considered in this manuscript, and an 
additional geometrical transformation (an extension of the original graph) has been introduced in Sections V C|V P" 



In this Appendix we demonstrate that the equivalence of the ice model and the Ising model to respective WBG 
models can be established without extending the original graph. To that end we extend the gauge group (i.e., the group 
of gauge transformations) by relaxing the constraints given by Eq. ( 15 1. Specifically, we relax the constraint iTat — i^ba 
in the definition of the edge variables, so that the edge variables tt^ are represented by pairs "Kf^^.b] = {^^ab, T^ba) of 
binary variables. In addition to the vertex weights fai'^a) we further introduce the edge weights gaiT^a)- We will also 
use a notation Qabi^^ab, i^ba) with a natural constraint gabi^^ab, T^ba) = gbai^^ba, T^ab) that makes this notation consistent 
with the original one. The partition function of such an Extended EBG, thus EEBG, becomes 



ZeEBG = X! n /"(TTa) W ga{T^a) 



(Bl) 



A gauge transformation of an extended model is described by a set {GabiT^abTT^'ab)} {a bjeSi local invertible 2x2 



matrices with no other restrictions. A gauge transformation for the vertex functions is still given by Eq. (14), whereas 
the edge functions are transformed according to the rule 



5{a,h}(7r{a,6}) g{a,b}i'^{aM) ^ X! '^"'^ ^^Cib.T^'ab)^^'' i'^ba,T^'ba)9abWab^'^'ba) 



•{a,b} 



where G denote the matrices inverse to Gab, i-e. 



X! ^"''(^afc' T^ab)Gab{T^ab, Kb) = ^(Kb^ Kb) 



(B2) 



(B3) 



Obviously the partition function of the extended model is invariant with respect to the extended gauge transformations. 

We further note that the EBG model (16 1 is a particular case of the EEBG model (Bl| that corresponds to the 
choice of the edge functions in a form 



gab{'^ab,T^ba) — S{TTab,T^ba) 



(B4) 



and the EBG-model gauge transformations, i.e., the ones that satisfy the constraints of Eq. (15), can be viewed as 
the extended gauge transformation that preserve the edge functions of the form, given by Eq. ( B4 ) . Obviously, any 
extended model is equivalent to an EBG model. 

The ice model represents a class of graph-based models, hereafter referred to as arrow models on graphs (AG), 
where a configuration is given by a graph orientation, i.e, a set of arrows associated with the non-oriented edges, 
rather than binary variables residing on the edges. Any AG model can be viewed as a particular case of EEBG model 
by associating with a local edge configuration of an AG model on {a,b}, determined by an arrow a ^ h, the edge 
configuration tTq^ with -Kab = and Tr^a = 1. Obviously the obtained EEBG model corresponds to the choice of the 
edge functions represented by the same off-diagonal matrix with (?a{>(0, 0) = (?af)(l, 1) = and 5o&(0, 1) = 3a6(lj 0) = 1. 
A homogeneous gauge transformation of this form is 



1 

71 



1 



(B5) 
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transforms the edge functions ga of an AG model into respective ones described by Eq. ( B4 1 . This estabhshes an 
equivalence between an AG model and respective EBG model. Note that the edge function for an AG model is 
symmetric and, therefore, can be viewed as a quadratic form. Then the gauge transformation can be interpreted as 
the quadratic form diagonalization. If performed over real numbers the result has a diagonal form with the signature 
(1,-1). The imaginary units in the matrix elements in Eq. (B5 ) are responsible for transforming the obtained diagonal 
matrix into the form given by Eq. (B4|. 



Consider a subclass of AG models, hereafter referred to as even AG models, whose vertex functions stay invariant 
upon a change the direction of all arrows, associated with the vertex. (Note that the ice model belongs to this 
subclass.) It is easy to show that the described above homogeneous gauge transformation transforms an even AG 
model into an even EBG model. Also note that in this even case the vertex functions of the result ing e ven EBG model 
are real, despite of the presence of the imaginary unit in the gauge transformation given by Eq. (B5). The suggested 
gauge transformation identifies the subclass of AG models (including, e.g., the ice model) resulting in WBG models 
upon the gauge transformation given by Eq. (arrow-to-EBG), and, therefore are easy. 

The Ising model on an arbitrary graph Q can be also viewed as an EEBG model with the vertex and edge function 
given by 



1, Vfe, C e 5g{a) TTab 

otherwise 



0, 



9ab{T^ab, TTba) 



7, TTab = T^ba 



(B6) 
(B7) 



which differ from Eqs. (27 1 and (28) actually by notations only. The gauge transformation that diagonalizes the vertex 
functions g^, viewed as quadratic forms, turns the Ising model into an even EBG model that belongs to the WBG 
subclass, as outlined in subsection |VD[ 

Finally, we note that the approach discussed in this Appendix, based on the extended gauge transformation, can 
be viewed as an alternative formulation leading to the well-known high-temperature expansion of the Ising model, 
when the expansion is interpreted as a loop series with only Z2-cycles providing non-zero contributions. 



